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ABSTRACT 

Satellite -to- satellite tracking and orbit computation accuracy 
is being evaluated on the basis of data obtained from near 
Ear*h spacecraft via the geostationary ATS-6. The near Earth 
spacecraft involved are Apollo-Soyuz, GEOS-3 and NIMBUS-6. 
In addition ATS-6 is being tracked by a new scheme wherein a 
single ground transmitter interrogates several ground based 
transponders via ATS-6 to achieve the precision geostationary 
orbits essential in satellite-to-satellite orbit computation. 

Also one way Doppler data is being recorded aboard NIMBUS-6 
to determine the position of meteorological platforms. Ac- 
curacy assessments associated with the foregoing mission re- 
lated experiments are discussed. 
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1.0 INTRODUCTION 

Orbit determination accuracy of . orbiting spacecraft is primarily a func- 
tion of the: 

• knowledge of the Earth's gravity field, atmospheric drag, solar 
pressure, solar and lunar gravity effects, 

• tracking system performance and validity of atmospheric refraction . 
corrections, and 

• tracking geometry and station location accuracy. 

The tracking system makes measurements such as range, range rate, angles 
and direction cosines to a spacecraft relative to a given tracking station. This 
data, typically on magnetic tape, is preprocessed (i.e. edited, pre-smoothed, 
changed to metric units, combined with calibration constants and so on) prior 
to being used as input to an orbit computation computer program. The quantity 
of tracking data required to achieve a desired orbital accuracy is directly linked 
with the spacecraft-to-tracking station relative geometry and dynamics. The 
orbit computation program is a means of mathematically providing a "best fit" 
(usually in a least-squares sense) of the measui’ed tracking data to the physical 
laws of orbital mechanics. 

The span of data used for an orbit determination may be as short as a few 
minutes, such as during critical launch or injection maneuvers, or data sam- 
pling over a period of several days may be used during high eccentricity earth 
orbits and planetary trajectories. 

Assuming reasonable tracking geometry the accuracy of spacecraft position 
and velocity determination will be primarily limited by tracking system per- 
formance for any computation spanning the data collection interval. That is, if 
continuous tracking is provided from a set of well surveyed stations the com- 
putation is essentially one of geometry. On the other hand the accuracy of 
orbit prediction based on an initial spacecraft vector determination will be de- 
graded as a function of time in direct relation to the accuracy to which physical 
parameters are modeled. This modeling includes gravitational fields, atmos- 
pheric drag and refraction effects, solar pressure, station location determination 
and so on. The most critical of these modeling parameters in terms of orbit 
determination accuracy is the gravity field model which at present is generally 
expressed in terms of a spherical harmonic expansion. 
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The question to be answered is what does one mean by "accuracy" and how is 
it to be measured? For the tracking system the accuracy of measurement is 
usually determined by means of static collimation tower tests where known 
distances, angles, and zero Doppler conditions are predictable to the limit of 
physical survey, time delay measurement, and frequency measurement. Pre- 
dicting tracking system performance under dynamic conditions is usually ac- 
complished by observing the behavior of associated sub-systems under sim- 
ulated signal conditions. Experience has shown that the tracking system 
accuracy can be specified, implemented and verified with a high degree of confi- 
dence (ref. 1). On the other hand it is much more difficult to verify the accuracy 
of the complex modeling required to predict orbits. There is a continuing 
effort underway to improve our knowledge of such fundamental modeled param- 
eters as the Earth's mass, the universal gravitational constant, speed of light, 
drag coefficients and so on. For it is easily shown that assuming perfect 
modeling the propagated (i.e. predicted) orbit errors due to tracking system 
uncertainties can invariably meet all current operational requirements for 
navigation, telemetry acquisition and satellite based Earth platform location. 

In this regard it is well established that satellite position errors of tens of 
meters will propagate into kilometers after only one or two weeks of orbit 
production. Such errors can only be reduced by increasing the number of track- 
ing observation periods or improving the Earth's gravitational field modeling. 

This paper presents a number of orbit error assessments as derived from 
satellite-to-satellite tracking and orbit computation involving the geostationary 
ATS-6 and the near Earth GEOS-3, Apollo-Soyuz, and NIMBXJS-6. 


2.0 TRACKING SYSTEM PRINCIPLES 

All of NASA's range and range rate radio satellite tracking systems consist of 
narrowband phase modulated signals where several ranging tones are sequen- 
tially modulated onto the carrier to provide a degree of range measurement 
ambiguity resolution. Once the ambiguity is resolved only the highest frequency 
tone is transmitted to assure maximum resolution. The high resolution range 
tone extends from 20 kHz to 500kHz depending on the specific tracking system 
under consideration. The range observation is basically a propagation time 
delay measurement from ground station to satellite and back. In the case of 
the current NASA satellite-to-satellite tracking experiments (e.g. ground station 
to the geostationary ATS-6 to near Earth satellite) four distinct propagation 
links are involved. 

The range rate relative to a given tracking station is generally observed by 
counting cycles of carrier Doppler where each cycle of Doppler corresponds 
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to a half wavelength change In total path length at the operating frequency. All 
satellite-to-satellite tracking is currently conducted at a nominal carrier fre- 
quency of 2GHz and hence each cycle of two way Doppler corresponds to 7.5CM 
of path length change. Strictly speaking the ->asie measurement is one of range 
change and only becomes a range rate when averaged over the Doppler counting 
interval. 


The foregoing is depicted in simplified form in Figure 1 which is appropriate 
for the satellite-to-satellite tracking associated with the geostationary ATS-6 
tracking of GEOS-3, Apollo-Soyuz and NIMBUS-6. Satellite-to-satellite track- 
ing permits the computation of near Earth satellite orbits by means of a single 
tracking station. NASA is implementing an operational tracking system of this 
type in the 1980 time frame by means of the geostationary Tracking and Data 
Relay Satellite System (ref. 2) . 






2.1 TWO WAY S A T£ L LITE ~ TO - S A TE L LITE TRACKING 


In two way tracking a signal is transmitted from a well surveyed ground station 
to a transponder which frequency translates the signal for re-transmission 
either to another spacecraft or directly back t*. :'he ground station. This paper 
deals primarily with satellite-to- satellite tracking and hence the ATS-6/ 
NIMBUS-6 tracking system (ref. 3 & 4) will be used as an example. Except for 
minor details involving the near Earth spacecraft transponder tins discussion 
also applies to the ATS-6/GEOS-3 and ATS-G/Apollo-Soyuz tracking. 

The tracking of a near Earth satellite via a geostationary satellite might be 
accomplished in a number of different ways. For example, tracking signal 
generation and data demodulation might be performed directly at the synchronous 
satellite and sent by telemetry to a ground station. The advantage of such a 
scheme is that inter-satellite measurements (synchronous to low satellite) can 
be separated from the total path delay. The principal disadvantage is that the 
relative complexity of a total ground station must be placed in orbit. Another 
possibility which is easier to implement electronically is 'the "bent pipe" con- 
cept where tracking signals between ground station and near Earth satellite are 
relayed back and forth via the geostationary satellite. In this scheme all track- 
ing data demodulation, digitizing and recording is performed at the ground 
station. The disadvantages is that total path delay is combined into a single 
measurement which in turn, adds a degree of complexity to the orbit determin- 
ation program. It might be thought that the geostationary satellite motion is 
negligible over the observation interval, however this is not the case. The 
reason is that geostationary satellites, are generally not maintained at zero 
degrees inclination. Nominal values of Inclination maintained for current NASA 
geostationary spacecraft extend from 1.5° to 6°. Such synchronous orbits are 
apparently more stable than zero degree inclination orbits. As a result the 
slant range relative to an observing ground station typically undergoes a sinu- 
soidal variation of several hundred kilometers over a 24 hour period. This 
effect will also be reflected in the corresponding ground to geostationary satel- 
lite range-rate (Doppler) measurement as tens of meters per second. 

The "bent pipe" scheme has been implemented for all NASA satellite-to-satellite 
tracking and orbit computation to date. The geometry of such tracking is shown 
in Figure 2. This concept has also been used in a "trilateration tracking" 
scheme (ref. 5) as indicated by Figure 3 to pinpoint ATS-6 while stationed over 
the U.S.A. at 94° W Longitude. Recent experience indicates that the success of 
satellite -to-satellited orbit determination is to a large measure dependent on 
the accuracy of the a priori estimates used for the geostationary satellite. 
Conventional one or two station tracking of geostationary satellites, while per- 
fectly suitable for most meteorological image registration and data acquisition 
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Figure 2. Basic Tracking Geometry 


purposes, is not generally adequate for satellite-to~satellite applications. 
Trilateration tracking provides the geometry and tracking data resolution 
which results in geostationary satellite a priori estimates at the 100 meter 
or better level. 

The same scheme is being used to track ATS-6 at its current position of 35°E 
Longitude with interrogation from the Madrid site and transponders at Madrid, 
Ascension Island and Johannesburg. 

• SYSTEM DESCRIPTION 

ATS-6 is in an earth-synchronous orbit at 35,800km. 
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Figure 3. ATS-6 Trilateration Tracking 


NIMBUS-6 is in a sun-synchronous orbit at 1100km. Figure 2 shows this basic 
geometry and the location of a ground tracking station. 

The vital ATS elements for tracking Nimbus and relaying Nimbus experiment 
data to the ground are a communications transponder and a 2 -GHz nine-meter 
parabolic antenna. The transponder translates 6 GHz ground station signals 
into 2 GHz signals sent to Nimbus and also translates the 2 GHz signals from 
Nimbus into 4 GHz signals which are sent to the ground. Figure 4 shows these 
frequency links between Nimbus, ATS, and the ground. The antenna has a 
nominal beam width of 1.4 degrees and a gain of 36 dB. It can be electronically 
scanned ±5 degrees off bores ight, and has monopulse capability to track the 
Nimbus satellite. However, the primary ATS antenna pointing mode for the 
tracking experiment is for the ground station to program the ATS with computed 
pitch and roll commands based upon ATS and Nimbus ephemeris data. 

Nimbus observes ATS by means cf an up-looking 2 GHz antenna array which 
has a nominal gain of 15 dB corresponding to a 3 dB beam width of 25 degrees. 
The system consists of the gimbaled antenna assembly, the gimbal drive elec- 
tronics, the power amplifier, transponder, and the digital electronics. The 
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Figure 4. Tracking Signal Power Levels and Pointing 


gimbaled high-gain antenna is directed to ATS by programmed activation of an 
X/Y mechanical mount. The antenna is located on the top of the Nimbus 
spacecraft. 

The tracking data recorded at the ground station consists of range and Doppler 
measurements in terms of time delays. The range data is elapsed radio wave 
propagation time. The Doppler data, which is a function of range rate, is 
recorded as the time required to count a fixed number of cycles of two-way 
Doppler. 

Figure 5 illustrates the overall tracking data processing required for this 
experiment. 
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Figure 5. Data Processing lor Satellite -to-Satellite Orbit Computation 

The overall ranging measurement consists of a measurement of the total round 
trip signal delay and involves: 

• The interrogation station located at F. osman, North Carolina (or 
currently from che transportable site at Madrid, Spain) 

• The phase-leaked ATS translation transponder and S-Band parabolic 
antenna, and 

• The Nimbus crystal-controlled transponder used in conjunction with 
the programmed S-Band helical antenna array (phase-locked trans- 
ponder for GEOS-3 and Apollo) 

The rate of range change ("average range rate") is observed as a Doppler 
shift and necessarily involves the relative motions of the two spacecraft and 
the ground station. 

The ground station typically will transmit at a 2 kw (CW) level, although 
it is capable of transmitting up to 10 kw. The highest resolution range 
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tone is lOOKHz with lower ambiguity resolving tones used during acquisition. 

The tracking signal generated at the ground station and transmitted to ATS is 
used for the coherent Doppler and tone ranging measurements. The signal 
generation is indicated in simplified form in Figure 6 (ref. 6). The use of the 
pilot carrier at 6150 MHz permits a coherent lockup with the ATS transponder, 
while the tracking signal (6137. S5 MHz in the case indicated) can be varied over 
a wide range to permit coherent tracking of other spacecraft or ground located 
transponders without reacquiring ATS in the frequency domain. The coherent 
tracking signal translation by the ATS transponder is indicated in Figure 7. 

The Nimbus translation transponder shown, in Figure 8 is interrogated by ATS 
at 2062.85 MIIz. The four-element antenna subsystem is directed to ATS by 
programmed activation of a gimbitled X/Y mechanical mount. The pointing 
information is normally loaded into the Nimbus memory via VKF radio link 
when the spacecraft is in view of a Nimbus command site such as the Fairbanks, 
Alaska, STDN. The antenna can also be controlled by direct access and control 
via the ATS- 5 command relay to Nimbus (S-Band). 

With reference to Figure 8 the incoming Doppler-shifted ATS signal is trans- 
lated by a frequency derived from the Nimbus 37.550 MHz crystal oscillator. 
This same reference is multiplied up to S-Band (2253.0 MHz) where it serves 
as the carrier for the Nimbus-to-ATS eight- watt link. The translated ATS-to- 
Nimbus signal is phase modulated onto the 2253.0 MHz carrier at a nominal 
modulation index of 1.5 radians. This system is based on the Goddard Range 
and Range Rate concept, where a crystal-controlled relatively broadband (sev- 
eral hundred KHz) transponder is employed. This system is made equivalent 
to coherent (i.e., phase-locked transponder) operation by proper ground station 
processing of the transmitted carrier which is coherent with the onboard ref- 
erence oscillator. The advantage of using such a transponder is that no fre- 
quency swept acquisition is required by the interrogating signal. The frequency 
excursion at 2 GHz due to one-way Doppler often approaches ±50 KHz, since the 
linear speed of near-earth spacecraft is typically 8 lan per second. 

Doppler data is the most accurate form of tracking data available for purposes 
of orbit computation because one cycle of Doppler is recorded for every half 
wavelength the spacecraft moves radially relative to the interrogating station. 

The group delay of the Nimbus transponder has been carefully calibrated over 
a wide range of frequency and temperature. Measured group delay repeatability 
is within ±20 nanoseconds with a nominal delay of 2.6 microseconds, Systematic 
one-way ranging errors introduced by the Nimbus transponder are thus expected 
to be less than ±3 meters. 
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Figure 6. Ground Station Signal Generation to ATS-6 
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Figure 7. Signal Translation by the ATS-6 Transponder 
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Figure 8. N. T MBUS-6 Transponder 
• TRACKING MEASUREMENT INTERPRETATION 

The "range" measurement is performed by comparing transmitted and received 
tone zero crossings, the highest resolution tone frequency in this case being 
lOOKHz. The "range rate" measurement is performed by counting a predeter- 
mined constant number of Doppler cycles and recording the time required to 
receive these cycles. Thus, both "range" and "range-rate" are recorded in 
terms of elapsed time. The raw data actually lists elapsed cycles of a 100 MHz 
clock; consequently the time readout is quantized to 10 nanoseconds. 

The highest resolution ranging tone used in this experiment is 100 KHz. Lower 
frequency tones are sequentially used during acquisition for ambiguity resolution. 
The lower tones are at 20 KHz, 4 KHz, 800 Hz, 160 Hz, 32 Hz, and 8 Hz. 

The tone ranging measurement is quite straightforward. However, its accuracy 
depends chiefly on the quality of preflight calibration of both the ATS and Nimbus 
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transponder group delay. Such preflight calibration data have been taken over 
a range of frequencies and temperatures. Indications are that with careful cali- 
bration the total systematic delay error in the ranging measurement can be 
held to a few meters of equ. . alent one-way range. 

The electronics for ATS- 6 satellite-to-satellite tracking have been so configured 
that the Doppler output is approximated by the following equation: 

- 2 f k _ _ _ 

f d = [a,fj + a 2 (rj + r 2 )] (1) 


where f d measured average Doppler frequency 
f t = uplink frequency = 6137.85 MHz 
k = 0.336 


aj and a 2 are scalar constants determined by equipment frequency 
multiplications 

r x = average range-rate ATS to ground site 

r 2 = average range-rate ATS to Nimbus, Apollo, or GEOS-3 


A detailed discussion of Doppler factors in satellite-to-satellite tracking is 
given in ref. 7. A description of the observations, data formats and system 
parameters associated with NASA-GSFC satellite-to-satellite tracking is given 
in ref. 8. 


In order to permit range; rate direction determination a fixed bias frequency 
(500 KHz) is added to the observed Doppler at the Doppler extractor. Thus the 
system counts the time, T c , required to accumulate N cycles of Doppler, f d , 
plus bias, f b . That is: 


T e = t t — (seconds) 

f b + f d 


where N is given in Table 1. 


(2) 
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Table 1 

NIMBUS Doppler Cycle Count 


N 

Data Sample Rate 

31995 

8 per second 

63990 

4 per second 

127980 

2 per second 

255960 

1 per second 

2559600 

6 per minute 


The GEOS-3 tracking uses a continuous Doppler count such that an accumulation 
or Doppler cycles over a count time T is given by: 

N 2-Nl=T(f b i-100f d ) {3 ) 


and the bias frequency is 2 x lo 7 Hz. For GEOS-3 the Doppler equation in the 
form of equation (1) is given by: 


f d = 2.247 x 10 9 




1.700\ 

2.247/ 



( 4 ) 


Here the scalars 2.247 and 1.700 are exact and c is the speed of light. 


ft 


• ATMOSPHERIC EFFECTS 

The range and Doppler measurements will also be biased by the Earth's tropo- 
sphere and ionosphere . Measurement biases of meters in range and tens of 
cm/sec in range rate can be expected at 2GHz. Atmospheric refraction effects 
can to a large extent be modeled out. Some of the work done in this area at 
NASA-GSFC is indicated in references 9, 10, and 11. 

The atmospheric range bias is frequency independent through the troposphere 
and inversely proportional to frequency squared through the ionosphere. The 
range rate bias, in addition to the foregoing, is proportional to the rate of scan 
through the atmosphere as well as to the magnitude of horizontal gradients . 
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• ORBIT DETERMINATION 


In the actual orbit determination program (NASA-GSFC Navigation Analysis 
Program) the four radio propagation paths are considered separately and the 
final orbit solution arrived at iteratively. This analysis piogram is a gener- 
alized least squares parameter estimation program designed to accept and 
process numerous types of tracking observations, and includes algorithms for 
rigorous treatment of single or multi-satellite time delay and delay rate meas- 
urements. The geopotential model can be selected from any of a number of 
available gravity field models. One such field currently used in this program 
is the Goddard Earth Model-2 (GEM-2) which is given in terms of spherical 
harmonics up to order and degree 22. Lunar and solar perturbations are pro- 
vided by means of the JPL ephemeris residing on permanent disk file at the 
NASA-GSFC IBM 360/95 computer facility. 


Studies have shown (ref, 12 & 13) that Earth gravit atio nal anomalies should be 
observable in near Earth satellite -to-satellite Doppler data. Preliminary 
analysis of ATS-6/Apollo-Soyuz data bears out the conclusions of these early 
studies ;ref. 14). 


2.2 NIMBUS-6 ONE WAY DOPPLER TRACKING 

One way Doppler measurements are being made aboard NIMBUS-6 in conjunc- 
tion with the satellite balloon and buoy tracking and meteorological experiment 
entitled "The Tropical Wind Energy Conversion and Reference Level Experi- 
ment (TWERLE). On the daylight portion of each orbit, when the NIMBUS is 
within range of the meteorological balloons (several hundred up at present) the 
Random Access Measurement System (RAMS), which is the space-borne segment 
of TWERLE, detects, demodulates and stores the one way Doppler shifted signal 
and sensor data transmitted by each platform (ref. 15). The power level of the 
beacon is on the order of 600 milliwatts at a nominal carrier frequency, f 0 , of 
401.2MHz. The expected uncertainty in this frequency is ± 5kHz and the maxi- 
mum one way Doppler shift at NIMBUS as calculated from -r/c f 0 is ± 10kHz. 
The onboard double conversion receiver derives successive translation fre- 
quencies of 345.6MHz and 55.575MHz from the NIMBUS 1.6MHz clock and con- 
sequently operates at a nominal intermediate frequency (IF) of 25kHz. The IF 
operating range is from 10 to 40kHz (ref. 16). Figure 9 indicates the principle 
of operation of the RAMS system. The NIMBUS operational orbit is computed 
by the NASA tracking network and the ephemeris (X, Y, Z and X, Y, Z) is used 
as computational input at the ground based computing center. Frequency offset 
(i.e. departure from nominal 401.2MHz Doppler) is solved for along with the 
beacon position coordinates Sj S n S 3 . Position can be determined in a single 
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R 2 = (X-Si) 2 + (Y-S 2 ) 2 + (Z-S 3) 2 

- dR • dR * dR* < x_s i) • ( y " s 2 ) • (2 _ S 3 ) . 

R = — X + — — Y + --Z = — X + -- Y + — Z 

dX <3 Y dZ R R R 

DOPPLER SHIFT AT SATELLITE = - -^-f T 

NOTE: GEOCENTRIC CARTESIAN 
EARTH FIXED COORDINATES. 
C= 2.997925 x 10 8 METERS SEC 

Figure 9. One Way Doppler Location of Beacon 

satellite pass. Balloon velocity is estimated from two successive passes 
(ref. 17). The software developed for satellite-to-satellite orbit computation 
(i.e. Navigation Analysis Program) has also been applied to the solution of the 
one way beacon problem and consistently recovers such parameters as fre- 
quency offset, timing bias and beacon location. The fixed reference platform 
used in experimenting with this data type (apart from the normal TWEELE 
Balloon operations) is located at the NASA Fairbanks, Alaska tracking station. 


3.0 EXPERIMENTAL RESULTS 

The following presents some recent results in the areas of satellite-to-satellite 
tracking, geostationary satellite trilateration and one way Doppler Earth fixed 
beacon location from a near Earth satellite. 
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3.1 SA TE LLITE -TO -S A TE LLITE ORBITS 


The geostationary ATS-6 spacecraft launched on 30 May 1974 has been the relay 
satellite for all NASA salellite-to-satellite tracking to date involving GEOS-3 
(launched 9 April 1975), Apollo-Soyuz (15 July - 24 July 1975), and NIMBUS-6 
(launched 12 June 1975). 

The expected error for the NASA range and range rate satellite-to-satellite 
tracking system is a function of many controlled parameters such as range tone 
frequency, sample rate, bandwidth settings, signal-to-noise spectral density 
ratios, spacecraft dynamics and so on (ref. 6). However the system is generally 
used with what might be termed a standard set of options such as; 100kHz maxi- 
mum range tone frequency, signal levels such that system is not thermal noise 
limited, 1 per second or 6 per minute data rate, and a 25Hz range tracking loop 
two-sided noise bandwidth. Table 2 lists the theoretical system performance 
for the foregoing selected options. Doppler averaging time is approximately 
one half the sample time interval for NIMBUS tracking and equal to the sample 
interval for Apollo and GEOS tracking. 

For averaging times, T, up to about 10 seconds the noise decreases as l/T. 

The principal Doppler noise contribution comes from receiver voltage controlled 
crystal oscillators and the analog to digital conversion. For longer integration 
times the Doppler noise is also influenced by noise falling off as 1/ (Fig. 10). 
This effect is attributed to the phase jitter in the transmitter reference signal 
used at the Doppler extractor. It should be pointed out that the least significant 
range bit recorded is 1.5 meters which is consistent with the best expected one 
way performance of 1.7 meters resolution. 

Table 2 indicates the predicted satellite-to-satellite tracking system measure- 
ment resolution. 


Table 2 

Tracking System Measurement Resolution 


Range (Meters) 

Range Rate (Cm/Sec) 

Systematic 

Random 

Systematic 

Random 

1.2 

• 

1.2 

Negligible 

0.03 



NOTE: 10 Sec. Averaging 

i 
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Figure 10. Range Rate Resolution 

Measured results indicate close agreement with expected system performance. 
System random errors or "noise" are generally observed by the least squares 
fitting of short data spans (i.e. 1 to 10 minutes) with polynomials of at least 5th 
degree to account for spacecraft dynamics. Care must be taken such that the 
polynomial itself does not introduce apparent error, (ref. 18). If the data is 
from a static or collimation tower test a least squares straight line fit is 
appropriate. 


Once an orbit is computed one indication of valid results is the difference be- I 

tween observed and calculated parameters over the data span. Such differences \ 

are often referred to as "orbit residuals". When the data is compared with the 

calculated orbit over one or more revolutions the results will be Indicative of 

uncertainty in modeled parameters such as imperfect gravitational field har- > 1 

monies and so on. It is clear that in order to improve such modeling through 

orbit solutions the basic orbit determination uncertainty introduced by tracking 

system errors must be less than the perturbation being solved for. Another j 
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test of validity is the independent determination of two or more orbits of a given 
spacecraft for the same time span epoch and differencing the position and veloc- 
ity components. Differencing independently determined oi’bits for a given space- 
craft over the same time span is often referred to as "orbit overlap analysis" 

(ref. 19). 

If independent tracking data sets are used to overlap the same orbit interval 
over which the overlap comparison is made the differences generally reflect 
tracking system performance, tracking geometry, data quantity, station location 
uncertainty and computational accuracy. However if one orbit is computed from 
a given data set and then predicted or "propagated" several days or more to 
overlap an orbit computed with a new data set, the orbit errors will be primarily 
due to uncertainties in gravitational field modeling. Orbit determination errors 
can thus be assessed either near the time of data observation or after orbit 
propagation. 

Finally, another measure of accuracy is the closeness of recovery of the same 
parameter or orbit using two or more independent means of comparison such 
as; ground tracking by lasers versus satellite-to-satellite tracking, or using 
different orbit computation programs to recover the same parameter and so on. 
Again the differences in results provide some measure of accuracy. One must 
be careful in interpreting such results since such comparisons usually imply 
one tracking system and/or computer program can be referred to as a "standard". 


• ATS-6/GEOS-3 ORBIT DETERMINATION 

Satellite-to-satellite tracking and trajectory computation analysis has been 
underway for approximately 6 months. ATS-6/GEOS-3 tracking provided the 
first example of this new data type. Preliminary indications are that the ac- 
curacy of the a priori position used for the geostationary satellite (ATS-6 in 
this case) is a very critical factor in achieving orbit solution convergence. 

The most effective procedure for two-way satellite-to-satellite orbit deter- 
mination appears to be as follows; 

(A) acquire satellite-to-satellite range sum and range sum rate data over 
several successive low satellite passes 

(B) obtain a reasonably accurate (i.e. position within several hundred meters) 
geostationary orbit by means of, for example, trilateration tracking 
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(C) obtain an approximate near Earth satellite a priori vector based on 
operational predictions 

(D) solve for both geostationary and near Earth satellites simultaneously. 

In this manner a 34 hour arc starting May 2 - 2300 Hours UT and ending 
May 4 - 0900 Hours UT consisting often passes of range sum and range rate sum 
satellite-to-satellite tracking data was processed to estimate the GEOS-3 orbit, 
the ATS-6 orbit, Targe data bias and solar radiation coefficients for each satel- 
lite. A priori estimates of each satellite epoch states were obtained from ground 
based tracking. Data rates were 6 per minute. The GEM-1 geopotential field 
was selected for this particular solution since this field has been most widely 
distributed to the user community. Orbit propagation studies are also currently 
being pursued using ATS-6/GEOS-3 satellite-to-satellite data in conjunction 
with evaluation of the GEM-1 through GEM-8 geopotential fields. The a poster- 
iori residuals of the fit were 1.3mm/sec R.M.S. for range rate sum residuals 
and 16 meters for range sum residuals (Fig. 11 & 12), The optimal weights for 
a priori estimates of states have not yet been determined. A solution when the 
ATS-6 satellite state is totally constrained yields larger residuals than the 
solution above. An unconstrained solution fails to converge. Preliminary over- 
lap tests indicate that satellite-to-satellite tracking can produce orbits com- 
parable to that derived from ground based tracking (ref. 19). For example the 
34 hour arc just described was used to produce two GEOS-3 orbits which were 
then overlapped. One orbit was computed from the first 24 hour data span and 
the second based on the last 24 hour span 2nd the overlap consisted of 12 hours. 
Comparison of the two GEOS-3 orbit vectors thus generated resulted in a mini- 
mum position difference of 10 meters and a maximum difference of 30 meters. 


These orbital differences represent the very first result of overlap as applied 
to satellite-to-sacellite tracking. Each of the two overlapped orbits was derived 
from data acquired during only three GEOS-3 revolutions. 


• ATS-6/APOLLO-SOYUZ OP BIT DETERMINATION 

The possibility of using high resolution satellite-to-satellite Doppler data in a 
short arc orbit computation tc detect gravitational anomalies led to the Apollo- 
Soyuz Geodynamics Experiment (ref. 14). The purpose of this experiment was 
to demonstrate the feasibility of detecting and recovering high frequency com- 
ponents of the Earth's gravity field by observing Doppler data from the near 
Earth (200KM altitude) Apollo satellite via the geostationary ATS-6 satellite. 
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Detectability haa been demonstrated through an analysis of the residual patterns 
in the satellite -to- satellite tracking data and comparing these patterns with pre- 
viously predicted signatures due to perturbations of gravity anomalies. The 
recoverability objective involves the actual estimation of the magnitude of these 
anomalies. 

The prime area of experiment data collection was the Indian Ocean Depression 
centered at 5° N. Latitude and 75° E Longitude. The experiment data collection 
phase was very successful. All data were obtained for the originally requested 
28 experiment revolutions. In addition, data were collected on 79 unscheduled 
revolutions. 

Preliminary results show that the detectability objective of the experiment has 
been demonstrated in the Indian Ocean Depression (Figure 13) area as well as 
over several other anomalous areas. Further analysis is required to demon- 
strate the possibility of actual recovery or estimation of the magnitude and dis- 
tribution of the anomalies. 

It should be pointed out that because of the short arc - typically 35 minute - 
orbit solutions used to detect anomalies, the results are essentially "independ- 
ent" of the particular spherical harmonic gravity field used. The GEM-1 field 
was used during the analysis which produced Figure 13. 


APOLLO CROSSING 
INDIAN OCEAN ANOMALY 
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3.2 ATS-6 TRIL A TER A T £ON ORBITS 

Experience to date indicates that well converged satellite-to-satellite orbit 
determination requires an a priori geostationary satellite position to an ac- 
curacy of a few hundred meters. This might be achieved by means of three 
or more widely spaced tracking stations. However a much more efficient 
means for accurate geostationary orbit computation has been established where- 
in a single^ station sequentially interrogates a number of widely deployed trans- 
ponder via the geostationary' satellite. The same station then records the track- 
ing data in exactly the same manner as during satellite-to-satellite tracking. 
Weekly trilateration tracking of ATS-6 is now being prrformed in support of 
the satellite-to-satellite tracking involving GEOS-3 and NIMBUS-6. The follow- 
ing ATS-6 trilateration tracking test was conducted during the checkout phase 
of the satellite-to-satellite tracking system. 


Geostationary Satellite Trilateration 


On 4 November 1974 a 24 hour ATS-6 trilateration tracking test (Fig. 3) was 
run using the tracking stations at Rosman N. C. and Mojave, California. The 
availability of essentially 2 separate data sets (Mojave and Rosman) covering 
the same 24 hour period made this test extremely valuable in assessing orbit 
computation accuracy since the prime transmit -receive sites (i.e. Rosman & 
Mojave) could be expected to contribute the major tracking measurement un- 
certainty. TIig transponders were located at Rosman, Mojave, NTTF Greenbelt, 
Maryland and Santiago, Chile. Each data stretch was approximatedly 5 minutes 
long and the data rate was one sample per 10 seconds. Over the 24 hour period 
each station sequentially tracked the transponders via ATS-6. Rosman and 
Mojave tracked over alternate 2 hour period, 0 . The data noise was at the system 
resolution level (Figure 14). The two position vectors calculated from the 
Mojave and Rosman data sets were overlapped and were in agreement to within 
20 meters at the center of the data span (Figure 15). It is reasonable that the 
center of the span should represent maximum accuracy since this point in time 
is bracketed by equal quantities of data providing maximum information regard- 
ing higher order time derivatives of position. As a matter of interest the 
effectiveness of the ATS-6 solar pressure modeling is indicated in Figure 16. 


3,3 NIMBUS-6 DOPPLER BIAS RECOVERY 

A one way Doppler satellite-to-satellite system clearly would be much simpler 
to implement that the two way' system described in this paper. The primary 
disadvantage of the one way Doppler sy'stem is that the near Earth satellite 
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1 .' ’ure 14. Tracking System Noise 




transmission frequency must somehow be accurately determined without direct 
measurement. Once the offset from nominal frequency has been determined 
one way Doppler satellite-to-satellite orbit solutions are possible. In all one 
way Doppler tracking systems used for position determination it is necessary 
to determine this offset between the nominal and actual transmitting frequency. 
The one way Doppler as recorded aboard NIMBUS-6 at 401.2 MHz has per- 
mitted an assessment of how well this frequency offset can be recovered in a 
single IS minute data pass. 

The NIMBUS-6 one way Doppler beacon location capability is being investigated 
in terms of Search and Rescue Applications. Single satellite passes over the 
Fairbanks, Alaska station are being analyzed. One question of interest is the 
recoverability of the frequency bias introduced by uncertainty in the ground 
beacon and spacecraft translation oscillators. Figure 17 indicates a typical 
Doppler pass. Figure 18 suggests that a frequency bias is quite observable in 
that one test using four distinct methods of calculation resulted in agreement 
to within 0.4Hz for a system with an expected measurement accuracy of ± 1Hz 
(ref. 16). The NIMBUS project calculation uses the algorithm described in 
ref. 17. The GTDS program is the principal NASA-GSFC operational orbit 
computation program. The NAP program is used primarily for orbital analysis. 
The geometric solution consisted of programming the observational equation 
indicated in Figure 9 and solving for the frequency offset along with the unknown 
position S x S 2 S 3 . All four cases used the same unadjusted NIMBUS-6 ephemeris 
which had been computed independently. 
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NOTE: FAIRBANKS STATION TO NIMBUS-6 
t 0 = 401.2 MHz 

START TIME 10 AUGUST 1975 
21 H 12 M 31.44 S 


Figure 17. Typical One Way Doppler Measurement 


• FAIRBANKS ALASKA BEACON 

• NOMINAL 401 .2 MHz 

• NIMBUS-6 RECEIVER 

• SINGLE PASS RECOVERY 

NIMBUS PROJECT 3098.5 Hz 
GTDS PROGRAM 3098.7 
NAP PROGRAM 3098.7 

GEOMETRIC 3098.3 

Figure 18. One Way Doppler Bias Recovery 
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4.0 CONCLUSIONS 

Satellite-to-satellite tracking system measurement resolution is predictable 
and readily verified. The resolution of the system described in this paper is 

0.05 cm/sec in range rate and 2 meters in range. 

A long arc (34 hour) GEOS-3 satellite-to-satellite derived orbit displayed RMS 
observed minus calculated values of 0.13 em/sec in range rate and 16 meters 
in range. Over the same arc the GEOS-3 orbit position differences as obtained 
from a 12 hour 3 pass overlap computation are 30 meters or less. 

The Apollo-Soyuz Geodynamics Experiment demonstrated that short arc satellite- 
to-satellite orbit solutions readily detect gravity anomalies. The Indian Ocean 
anomaly Doppler signature (observed minus orbit calculation) was at the expected 
level of 5 cm/sec. 

Experience to date indicates that successful satellite-to-satellite orbit deter- 
mination is to a large extent dependent on the accuracy of the a priori geo- 
stationary orbit. Trilateration geostationary satellite tracking to consistently 
reduce ATS-6 position errors to less than a few hundred meters is being per- 
formed regularity in support of all AT8-6/GEOS-3 and ATS-6/NIMBUS-6 
satellite-to-satellite tracking currently being performed. 
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